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We consider the higher-order Fierz transformation, which corresponds to expanding a product 
of t/jFip terms into a sum of products of Dirac densities and currents. It is shown that the Fierz 
transformation can be obtained by solving a large system of linear equations with fractional complex 
coefficients, which is practical at least up to fourth power. 
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I. THE PROBLEM 



Fierz transformation [[l] is a name given to the expression of a certain product of nondiagonal matrix elements of 
Dirac T-matrices as an expansion into products of diagonal matrix elements, such as 
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(ar^)(El» = c u {aT k a)(J)Tib). (1) 

k,l=l 

> , 

Here stands for one of the 16 Dirac matrices {1, 75, 7^, 757^, cr^} constituting a linearly independent basis in the 
space of complex 4x4 matrices. The matrix elements denote products in Dirac-index space only, i. e., 

r-> 
o 

[ The Fierz transformation is useful for expressing exchange matrix elements in terms of densities, currents, and other 
diagonal ones, which greatly eases their use in, for example, relativistic mean-field theories. For this reason it has 
\ been studied extensively, see, e. g., the papers by Y. Takahashi || or generalizations to SU(n) (J]. 

In the context of nonlinear self-coupling of meson fields, higher-order versions of the Fierz transformation have 
become of interest (for a review see |b)). If we express the order as the number of T-matrices involved, the above Eq. 
(El) is of second order, and in this paper we will be concerned with constructing expansions in third 



(ar i 6)=^ a (r,t)r i ^ 6 (r,t). (2) 
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and fourth order: 
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{aT i b){bT j c){cT k a) = ]T (aTia)(bT m b)(cT n c) (3) 

l,ra,n— 1 
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(aT l b)(br j c)(cT k d)(dTia) = c mnpq (aT m a)(bT n b)(cT p c)(dT q d). (4) 

m,n,p,q— 1 

Here, a general computer-based solution method will be discussed, using cither computer algebra systems such as 
Mathematica j|] or conventional programming allowing exact manipulation of fractional complex numbers. Although 
this approach is quite general, it will be applied preferentially to the case of symmetrized matrix elements with only 
the matrix 1 appearing in the expression to be expanded. 



II. SYMMETRIZATION 



In many applications, the wave function indices will all be summed over, so that it is sufficient to deal with an 
expression symmetrized over the indices. Thus, on the left-hand side of Eq. (|l|) in the second-order case we can write 

(a r, b) (b r 3 a) = \ ^ [(a r, b) (b r s a) + (br ia ) (a r, b)] , (5) 

ab ab 
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and the right-hand side of Eq. (|l|) will be symmetrized exactly in the same way. 

Using the notation X){a&c } ^° re f er to the sum over all permutations of the symbols a, b, c . . ., we can reformulate 
the Fierz transformation problem for the symmetrized matrix element as 

{ah} {ab} kl 

where the factor I has been dropped on both sides. 

It is important to realize that because of the product structure, the right-hand side is symmetric under an exchange 
of the r matrices as well. It is thus useful to introduce a notation for symmetrized terms, 

r fc ®r z = (ar fc a)(6r z 6) = {aT k a)(bTib) + (aT ia )(bT k b), (7) 

{ab} 

which can easily be generalized to higher order, for example, in third order the symmetrized problem becomes 

^ (aib)(bj c)(cka) = ^ c imn^i® (8) 

{abc} l<m<n 

with 

r, ® r m ® r n = ^ (ar,a)(&r TO &)(gr n c), (9) 

{abc} 

where symmetrization could equivalently be carried out in the indices l,m,n instead of a,b,c. Terms of fourth and 
higher orders are defined analogously. 

III. ALGEBRAIC DETERMINATION OF THE EXPANSION 

The second-order Fierz transformation as defined in Eq. (|]J) can be viewed as a system of equations obtained by 
comparing coefficients in the 4 4 = 16 2 dimensional space spanned by the spinors ip a , ipb, "0a, and ipt,. The coefficients 
are given by the components of the T-matrices and thus can be expressed as complex integers. The unknowns Cki are 
16 x 16 in number, so that we have exactly the right number of equations, and since the T-matrices form a basis for 
the 4x4 complex matrices, the decomposition ([!]) is always possible. 

The solution of this system of linear equations can be carried out using the standard Gauss elimination algorithm, 
provided that the coefficients are not treated in floating point arithmetic, but as exact complex fractions. 

For third-order Fierz transformations the dimension of the system of equations is 16 3 and it is 16 for fourth order, 
that is, the complexity in going from second order to fourth order is in the ratio 1:16:256. For the latter case practical 
solution would require substantial computing resources, but fortunately in cases of practical interest the number of 
terms in the expansion can be reduced substantially by symmetry and invariance requirements. In the symmetrized 
case of the preceeding section, for example, the dimension in fourth order is reduced by the number of permutations 
4!. 



IV. SELECTION OF THE TERMS IN THE EXPANSION 

The expansion into products of the diagonal matrix elements of the T-matrices is always possible, but usually is 
not the most useful expression of the Fierz transformation. To see this, let us look at an important special case: that 
of identity matrices on the left-hand side. The decomposition problem thus is 

(ab)(ba) = J2 c M^j a)(bT k b), (10) 

jk 

or, in symmetrized form, 

^(ab)(ba)=Y / c jk T j ®T k . (11) 

{ab} jk 
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Since the left-hand side is a Dirac scalar, this means that the right-hand side also can contain only scalar combina- 
tions of T-matrices. The only scalar combinations built out of products of two T-matrices are 1 <g> 1, 75 <§) 75, 7^ (8 7^, 
757m <8 757 M , an d cr^y ® a^ v , assuming the familiar index summation convention. The Fierz transformation problem 
in this case thus can be restated as (note that here because of the complete symmetry of all terms, the symmetrization 
can be omitted): 

(a b) (b a) = a (a a) (6 b) + c 2 (a 75 a) (6 75 b) + c 3 (a 7m a) (6 7^ b) (12) 
+04(0757,, a) (6 757^ b) + c^aa^ a)(ba' 11 ' b). 

Note that symmetrization works slightly differently in this case: the scalar products sometimes make certain index 
combinations appear repeatedly in the expansion, but it is still sufficient to include only one ordering of the T matrices 
in the symmetrized terms. 

Eq. (Jl^) corresponds to 256 equations for the 5 unknown coefficients. Clearly most equations will be redundant; 
eliminating them from the statement of the problem, however, turns out to complicate the solution, but the high 
degree of redundancy provides a welcome check for completeness and consistency of the assumed decomposition. 



V. SOLUTION FOR THE FOURTH-ORDER CASE 

Equation (|l^) will be used to illustrate the method of solution. Actually two different approaches were used 
depending on the complexity of the problem. 

For second and third order a very simple but flexible approach was implemented in Mathematica. The spinors were 
expressed as four-component vectors containing symbols of the form 

ip a ->• (al,a2,a3,a£), ip a -> (aal, aal, aa3, aaA), ip b -> (61, 62, 63, 64), ^-^(661,662,663,664), (13) 

The scalar products with the T-matrices can then be evaluated straightforwardly, yielding a representation of Eq. ( |T^ ) 
as a linear equation in the coefficients Cj with coefficients biquadratic in the spinor components. The symmetrization 
is done automatically. The coefficients are determined successively by choosing one term in which c\ occurs; if this is, 
e. g., ci aal a2 662 64, the coefficient of aal a2 662 64 is extracted from the equation, yielding a linear equation in the 
Ci alone, which is solved for c\. This is inserted into the equation, reducing the number of unknown coefficients by 
one, and the process is repeated until all ci have been found. If any terms are then left in the equation, the expansion 
was not complete. 

While this is a very straightforward and not particularly efficient solution, the steps can be automated using built-in 
functions, and it is quite flexible, since new terms can be added by simply writing them down as symbolic expressions 
in the spinors and T-matrices. For the fourth order, however, this process turned out to be too inefficient, so that as an 
alternative a Fortran-90 code was developed that uses a data type for fractional complex numbers and straightforward 
Gauss elimination. Programming the individual terms, though, requires substantially more coding. 

Specifically, the four independent components of each spinor are represented by a number 0, 1, 2, 3, corresponding 
to a two-bit integer. The four spinor indices are then combined into an 8-bit index, which indicates the index of 
the equation to which this combination contributes. For example, the expression (a 6) (6 a) has (as only one of the 
nonvanishing matrix elements) a value of one for the 3-component of spinors ip a and ipb and the 2-component of 
spinors ipb and ip a . If we arrange the index with ip a , ip a , ipb, ipb in order of decreasing magnitude, the index for the 
equation will be 

3x 4 3 + 2x 4 2 + 2x4 + 3= 11101011 2 = 235. (14) 

The program will loop through all spinor component combinations and add the generated coefficients at the corre- 
sponding position in the system of equations. 

To help with the proper selection of terms in the expansion, the code checks whether any of the coefficients of the 
term being generated are nonzero, since symmetry may cause cancellation in a non-obvious way, and also whether a 
term is directly proportional to one previously generated. 

The system of equations is then solved using Gauss elimination, where more general linear dependences will become 
apparent. If the elimination does not solve all of the equations, it is clear that the decomposition is not complete and 
this is then indicated. 
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VI. RESULTS 



As sample results, we give the decomposition of the symmetrized term in second, third, and fourth order for the 
case of identity matrices, as tables of terms with the corresponding expansion coefficients. For the higher-order terms, 
some remarks about new features are made. 

A. Second order 

The result for second order is 

(ab)(ba) = l(aa)(bb) + l(a l5 a)(b l5 b) + l(a^a)(brb) (15) 
- 2 (o 757^ a) (b 757^ b) + f (a a) (b a»» b) . 

This result is repeated in Tabic | in order to indicate the correct reading of the tables for higher order, where, however, 
a full symmetrization of both sides of the equation becomes necessary according to Eqs. (g) and 

B. Third order 

In third order the symmetrization is no longer trivial as it was in the second-order case. It may be surprising that 
terms with 7507,,, must be included; these can be equivalently formulated using the identity 

750>„ = ^ k \hu<^ kX , (16) 

but retaining the matrix 75 makes the space-reversal properties of the terms more readily apparent. Note that either 
way our basis still consists of only 16 linearly independent matrices. The resulting transformation is given in Table 

B 

C. Fourth order 

In fourth order the number of possible terms becomes larger and it is difficult to see which are independent. Using 
the basic building blocks 1, 75, 7^, 757^, er^, and 7507^ in all possible combinations fulfilling the condition of 
coupling to a Dirac scalar, which also implies an even number of 75-matrices, proves sufficient, but also leads to many 



dependent terms which are eliminated with the program's help. The final result is given in Table III 



VII. OTHER APPLICATIONS 

A similar method can be applied to nonsymmetric and nonscalar terms; in this case of course more terms will 
appear in the expansion. They can be constructed as before from the basic building blocks 1, 75, 7^, 757^, a^, and 
75C/J!/ in a way that yields the desired Lorentz transformation properties. 

As an example, the following decomposition was obtained: 

(a a^ v b) (b 7" a) = 

-\{aY a) (b b) - ± (a a) (b Y b) 

+^(a^a)(bb)-^(da)(b^b) (17) 
+f( s 757^ a) 75 b) + I (a 75 a) (6 75 7 M 6) 
-4(0 757^ a)(bj 5 a^ b) + 2(0750)^ a){b^ u b). 

Those terms where the exchange of indices a and b changes the sign will drop out in the symmetrized version of this 
result, which is 

\ [(a «v b) (b Y a) + (b a) (a 7" 6)] = -\ (a ^ a) (b b) - \{a a) (b 7" b) (18) 

+ T 757m «)(> 75^) + §(675 a) (6 757^6). 
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VIII. SUMMARY 



Once it is realized that constructing a Fierz transformation of any order essentially just means solving a linear 
system of equations with fractional complex coefficients, standard techniques of numerical analysis are sufficient to 
solve the problem. The calculation is now feasible up to fourth order, and the main work that remains is the expression 
of the transformation through meaningful couplings of the T matrices, which still requires treating each transformation 
separately. Our work here demonstrates that tractable solutions may be possible for each of these cases. 

Commented Mathematica and Fortran-90 programs are available from the authors. 
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TABLE I. Terms and expansion coefficients appearing in 
second order. 

Term Coefficient 

1®1 3 

75 ® 75 \ 

7m ® 7 M \ 
757m® 7s 7 M ~\ 



TABLE II. Terms and expansion coefficients appearing in 
third order. 

Term Coefficient 

101181 Ye 

1 <g> 75 <8> 75 je 

1 ® 7m ® 7" i| 

1®757m®757 m -j| 

1 <g> fT M „ ® ^ 

7m ® 7s7^ ® 75^ M " -f 

75 ® ® -y 5 a^ ^ 



•5 



TABLE III. Terms and expansion coefficients appearing in 
fourth order 



Term Coefficient 



1 



1 75 ® 0>f ® 

' In ® 757^ ® 750" Ml/ 

75 ® 75 ® 75 ® 75 

75 ® 75 ® 7m ® 7 M 

75 ® 7s ® 7s7m ® 757 M 



75 ® 75 ® M^ 



75 ® 7m ® 757^ ® a 



7m ® 7 ® 7w ® 7' 
7m ® 7 M ® 757^ ® 757" 
7m ® 7" ® 7s7m ® 757" 



J_ 

64 

_3_ 

32 



1 ® 75 ® 75 

1 ® 7m ® 7 M H 
1®7s7m®757 m -55 
1 ® a uv ® (7 Ml/ 



_3_ 

64 

M^ _3_ 

32 



_1_ 

64 

J_ 

32 

J_ 

"32 



uv _3_ 

' 64 



64 

__3_ 

32 

_1_ 

16 



7a ® 7 A ® ff/." ® el 

Ak _j_ 

16 



7a ® 7 ® "V*" 

757m ® 757 M ® 757^ ® 757 

7s7a ® 757 A ® f m^ ® c 

757a ® 757" ® ® crAl ' n 

3 

256 



64 

__3_ 

64 



64 
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